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Introduction 


This  paper  outlines  a  theoretical  study  of  the  electronic  structure 
of  several  transition  metal  organometal 1 ic  compounds  having  the  structure. 
(H3P)2MXY  (X,Y  m  H,  CHj)  carried  out  using  the  self-consistent-field-Xa  . 
scattered  wave  (SCF-Xa-SW)  method.  It  examines  the  characteristics  of  the 
upper  occupied  molecular  orbitals  of  these  complexes,  and  looks  for  correla¬ 
tions  between  these  characteristics  and  the  relative  rates  of  reductive 
elimination  of  the  elements  of  XY  from  them  (eq.  1,  L  ■  PH^).  These  reactions 
and  the  reactions 


L 

|  X  H  -  Pt,  Pd,  Ni 

I*  +  Y  X,Y  -  H,  CH3 

L 


(1) 


which  are  their  microscopic  reverse,  are  models  for  the  individual  steps  in  a 
number  of  useful  processes  catalyzed  by  transition  metals  (homogeneous 
hydrogenation,  coupling,  and  reduction;  perhaps  heterogeneous  hydrogenation 
and  dehydrogenation.)^ 

The  calculations  described  here  were  stimulated  by  two  qualitative 

general iza tons  which  are  emerging  from  current  experimental  work  in  the 

g 

organometal lie  chemistry  of  the  d  metals.  First,  the  rates  of  the  reduc¬ 
tive  elimination  reactions  of  the  type  represented  by  eq.  1  appear  to  de¬ 


crease  in  the  order  k^  >  ^(h)^  >  JL^ch^  for  a  series  in  which 
the  metal  and  its  ligands  l  are  constant.  Second,  the  rates  of  these  reactions 


seem  to  Increase  In  the  order  <  ^NKCH^  in 

2-4 

a  series  in  which  only  the  metal  changes.  Unfortunately,  there  are 

presently  few  strictly  comparable  examples  in  either  series,  and  the  value 

of  these  general izat ions  remains  to  be  established.  Further,  we  must  assume 

that  those  reductive  eliminations  which  are  known  follow  the  same  mechanism 

and  are  concerted  in  order  to  compare  their  rates. These  assumptions  are 

probably  correct  for  many  cases,  but  convincing  experimental  evidence  sup- 

2-4  6 

porting  their  correctness  exists  in  only  a  few.  *  Their  uncertain 
generality  not  withstanding,  these  two  series  serve  as  a  starting  point  for 
theoretical  consideration  of  this  type  of  chemistry. 

At  the  outset,  we  acknowledge  that  this  work  is  directed  to  only 
one  of  several  factors  which  might  in  principle  influence  the  relative  rates 
in  these  series:  viz. ,  the  electronic  structure  of  the  I^MXY  group.  Entropic 
effects  (which  are  probably  similar  for  all  the  reactions  subsumed  by 
eq  1^)  and  nonbonded  steric  effects  (which  may  be  significantly  different  for 

g 

different  of  these  reactions  )  are  neglected  entirely.  We  touch  only 
briefly  on  differences  in  the  energies  of  the  assumed  products  (L^M ,  X-Y)  , 
and  make  no  estimates  of  energies  of  transition  states.  Limitations  intrinsic 
to  the  version  of  the  SCF-Xa-SW  method  used  here  further  constrain  the  types 
of  questions  concerning  local  electronic  structure  which  can  be  addressed 
(see  below).  The  objective  of  the  work  is  thus  not  to  discuss  which  features 
of  the  electronic  structure  of  the  metal  and  Its  directly  bonded  ligands 
determine  the  relative  rates  of  eq  I,  or  even  to  try  to  establish  whether 
this  electronic  structure  does  in  fact  dominate  these  rates.  Addressing 
either  question  would  require  the  ability  to  construct  potential  surfaces 
(including  enthalpic  contributions  due  to  nonbonded  interactions  and  entropic 
terms)  describing  the  decompositions  in  a  detail  which  is  presently  beyond 
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the  capability  of  this  (and  other)  theoretical  methods.  Instead,  the 
objective  of  this  work  is  to  explore  the  limited  question  of  whether 
the  trends  in  rates  correlate  wi th  physically  interpretable  characteristics 
of  the  local  electronic  structure  of  the  reactants,  L2MXY.  If  such  corre¬ 
lations  exist  and  can  be  identified,  they  may  be  useful  in  guiding  the 
synthesis  of  mechanistically  or  synthetically  useful  organometal I ic 
complexes.  Any  more  sophisticated  effort  to  distinguish  between  "correlation 
and  "cause"  (n  this  area  of  organometal 1 ic  chemistry  must  be  postponed  until 
more  and  better  experimental  data  (especially  data  establishing  the 
structures  of  typical  transition  states)  and  more  accurate  theoretical 
methods  are  available. 

Previous  theoretical  discussions  of  rates  of  reductive  elimination 

Q  1  A 

reactions  include  works  by  Hoffmann,  et  al.  ,  Braterman  and  Cross  ,  and 

Akermark  .  Additional  relevant  calculations  have  been  described  by 

12 

Goddard,  Anderson,  and  others  . 

Methods 

Calculations  were  carried  out  using  the  SCF-Xa-SW  method  described 

elsewhere'^.  The  general  strengths  and  weaknesses  of  this  method  have  been 

14 

summarized  by  Salahub,  et_ aK  Here  we  emphasize  only  points  which  are 
relevant  to  the  particular  problem  being  explored.  The  SCF-Xa-SW  method 
has  two  strengths  vis  3  vis  other  theoretical  methods.  First,  it  can  cal¬ 
culate  wavefunctions  and  relative  energies  of  orbitals  for  many-electron 
systems  economically,  and  provides  good  wavefunctions  for  species  charac¬ 
terized  by  pronounced  ionic  character  in  bonds.  It  is  thus  practical  to 
use  this  method  to  examine  organometal 1 ic  compounds  containing  heavy 
elements  such  as  platinum,  palladium,  and  nickel.  The  ability  to  calculate 
electron  distributions  in  orbitals  is  especially  useful  in  identifying 
those  orbitals  which  contribute  significantly  to  the  bonding  (or  anti- 


bonding)  of  the  metal  to  its  ligands.  This  type  of  identification  is  not 
always  obvious,  especially  when,  as  in  the  complexes  examined  here,  bonding 
is  distributed  among  several  orbitals,  rather  than  localized  in  one.  Second, 
the  eigenvalues  generated  in  a  calculation  can  be  related  directly  to  orbital 
electronegativity  and  provide  a  characterization  of  the  electronic  structure 
which  is  obtained  with  more  di ff icul ty  using  other  methods.  The  reiat ive 
energies  of  orbitals  (that  is,  the  differences  in  energy  between  them)  are 
calculated  relatively  well  using  the  SCF-Xa-SW  method,  and  one  of  the  most 
useful  applications  of  the  method  has  been  in  the  calculation  of  electronic 
absorption  spectra  for  transition  metal  complexes  (using  the  "transition 
state"  method) 

The  present  SCF-Xa-SW  calculations  have  the  important  weakness 
that  they  give  unreliable  total  electronic  energies  of  molecules.  They  are 
thus  unsuitable  for  the  calculation  of  potential  surfaces,  and  for  problems 
which  require  energy  minimization.  The  type  of  calculation  employed  in  this 
paper  is  not,  for  example,  capable  of  estimating  ground-state  geometries^, 
and  certainly  cannot  be  used  to  estimate  structures  or  energies  of  transition 
states.  It  cannot,  therefore,  calculate  the  quantity  of  most  direct  interest 
in  any  study  of  reaction  rates  —  that  is,  the  difference  in  energy  between 
ground  and  transition  states  (the  activation  energy).  Efforts  to  use  the 
SCF-Xa-SW  method  at  this  level  to  examine  problems  in  kinetics  must  thus 
rely  on  inferences  based  on  the  (known  or  estimated)  ground  state  structures 
of  the  reactants  and  products,  and  on  the  energy,  symmetries,  and  shapes  of 
the  orbitals  of  these  materials.  Orbital  correlation  diagrams  have  been 
employed  as  an  aid  in  similar  analyses  employing  other  types  of  molecular 

9 

orbital  calculations.  As  we  shall  show,  and  In  agreement  with  previous  work  , 
they  are  relatively  uninformat Ive  in  the  particular  problem  chosen  here. 


In  order  Co  use  orbital  correlation  diagrams  in  any  fashion,  it  is 
necessary  to  be  able  to  calculate  orbital  energies  with  useful  accuracy. 

The  fact  that  the  version  of  the  SCF-Xa-SW  program  used  here  can  estimate 
individual  orbital  energies  usefully  but  that  it  cannot  estimate  total 
electronic  energies  deserves  some  explanation.  In  order  to  facilitate  rapid 
computation,  an  additional  approximation  is  introduced  in  calculating  the  Xa 
total  energy  that  is  not  made  in  calculating  individual  orbital  energies: 
the  charge  density  generated  from  the  one-electron  wavefunctions  is  spherical  1 
averaged  within  the  atomic  sphere  boundaries  and  assumed  to  be  a  constant  in 
the  region  between  the  spheres.  This  procedure  introduces  major  errors  in  the 
value  of  the  total  energy.  Dramatic  improvements  of  the  total  energy  have 
been  achieved  when  the  full,  unaveraged,  charge  density  generated  is  used 
in  evaluating  the  Xa  total  energy1^.  These  improvements  suggest  that  the 
discrepancies  between  total  energies  calculated  in  the  Xa-SW  model  and 
experiment  are  caused  mainly  by  the  approximation  of  a  spherical  charge 
distribution  and  are  not  a  reflection  of  intrinsic  defects  in  the  Xa  method. 
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These  distances  and  angles  are  unfortunately  poorly  based  in  detailed 

experimental  precedent,  but  are  sufficiently  accurate  to  examine  qualitative 

trends  in  the  complexes  involved^.  In  any  event,  these  choices  resemble 

q 

closely  those  used  by  Hoffmann,  et_  ak  in  a  useful  study  of  reductive 
elimination,  and  facilitate  comparison  between  the  results  of  these  studies. 
Caicul at  ions 

M A# 

1 8 

Schwarz's  values  were  used  for  the  atomic  exchange  parameters, 

19 

except  for  hydrogen,  for  which  0.77725  was  used  .  For  the  extramolecuiar 


and  intersphere  regions,  a  weighted  average  of  the  atomic  a's  was  employed, 
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the  weights  being  the  number  of  valence  electrons  In  neutral  atoms.  Over¬ 
lapping  sphere  radii  were  used^'^'. 

To  provide  a  check  on  the  consistency  of  the  parameters  used  for 
carbon  and  hydrogen,  the  transition  state  method  was  used  to  calculate  the 

first  ionization  potentials  for  dihydrogen,  methane  and  ethane.  These  values, 

22 

together  with  experimental  values  ,  are  (eV):-^,  17.6,  15.4;  CH^,  12.7, 

12.6;  C^H^,  11.0,  11.5.  A  similar  calculation  for  the  metals  gives  (calc'd., 

obsd. ) :  Pt,  8.77,  9.00;  ?d,  8.57,  8.33;  Mi,  7.89,  7-6 322,23.  Using  the  same 

method,  Parr,  et_  a]_.  have  calculated  the  electron  affinities  for  the  elements 

24 

up  to  atomic  number  54  with  a  similar  level  of  agreement  .  Comparison  of 
these  several  experimental  and  calculated  values  suggest  the  values  calcu¬ 
lated  by  the  SCF-Xa-SW  method  to  be  uncertain  by  about  1  eV,  We  place  this 
level  of  agreement  between  calculated  and  experimental  numbers  in  perspective 
by  noting  that  it  is  unlikely  that  the  total  spread  in  values  for  AG  for 
the  reductive  elimination  of  XY  from  the  LjMXY  complexes  studied  here  exceeds 
1  eV.  Thus,  the  uncertainty  in  quantitative  precision  of  the  computational 
method  is  approximately  the  same  as  the  magnitude  of  the  phenomenon  being 
studied;  accordingly,  only  qual i tative  trends  in  the  calculations  are 
interpretable. 

Atomic  charges  were  calculated  using  a  qualitative  procedure,  and 

25 

are  useful  only  for  comparisons:  the  absolute  values  are  net  reliable 
The  intersphere  and  extramolecular  charge  was  partitioned  among  the  atomic 
centers  in  proportion  to  the  intrasphere  charges.  That  is,  atomic  charges 
were  estimated  by  normalizing  the  total  number  of  valence  electrons  within 
all  the  atomic  spheres  to  the  total  number  for  the  molecule.  The  sphere 
sizes  for  a  particular  element  also  vary  slightly  from  compound  to  compound; 


1 


this  variation  also  complicates  quantitative  interpretation  of  calculated 
charge  distributions. 

Resul ts 


L2PtH2'  L2Pt^CH3’  and  L2Pt^CH3^2*  1  gives  calculated 


^  ^  ^  ^  ^  <v  ^  ^  ^  /w  /V  ^ 


eigenvalues.  Only  orbitals  involved  in  bonds  between  P,  Pt ,  and  X(Y) 

(and  nonbonding  d  orbitals  on  Pt)  are  shown;  inner  shell  orbitals,  and 
P-H  and  methyl  C-H  bonds  are  omitted.  Figures  2-4  summarize  the 
composition  of  Important  molecular  orbitals  ("important"  is  defined  by 
the  orbital  correlation  diagrams  as  orbitals  involved  in  bonds  which  are 
made  or  broken  during  elimination  of  XY  from  l^PtXY;  see  below). 

For  comparison,  Figure  5  gives  similar  information  for  L^Pt(O). 
Reductive  elimination  of  XY  from  l^PtXY  with  generation  of  L^PttO)  is 
certainly  accompanied  by  a  spreading  of  the  L-Pt-L  angle.  We  cannot 
estimate  the  value  of  this  angle  at  the  transition  state.  Figure  5  arbi¬ 
trarily  gives  orbitals  for  an  L-Pt-L  bond  angle  of  141°;  Figure  6 
indicates  the  change  in  these  orbital  energies  with  angle. 

Classification  of  molecular  orbitals  for  L^PtXY  as  bonding, 
nonbonding,  or  antibonding  by  Inspection  of  the  wavefunct ions  was  not 
always  straightforward,  for  two  reasons.  First,  certain  wavefunct ions  contain 
significant  contributions  from  platinum  6s  and  6p  atomic  orbitals,  and  the 
internal  nodes  of  these  orbitals  complicated  their  interpretation.  These 
internal  nodes  are  qualitatively  most  obvious  around  platinum  in  the  4^ 
orbital  of  (H^P^PtHj  (Figure  2)  and  are  due  to  the  6py  orbital.  Although 
there  are  obvious  nodes  between  Pt  and  both  H  and  P  in  the  electron  density 
plot  for  this  orbital,  the  Pt-H  and  Pt-P  interactions  are  bonding.  In  other 
orbitals,  it  was  not  always  possible  to  identify  the  origin  of  the  nodes 
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both  C  and  C  symmetries  revealed  Insignificant  differences.  The  figure 


Figure  5.  Orbitals  for  L2Pt(0)  for  angl 
Fl^'jre  2  contains  nomenclature. 
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Figgre  6.  Eigenvalues  for  L 
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by  Inspection.  Second,  for  orbitals  containing  a  significant  contribution 
from  the  platinum  5dz2  atomic  orbital,  It  Is  possible  to  have  coexisting 
bonding  and  antibonding  contributions  reflecting  overlap  of  the  ligand 
orbital  with  the  torus  or  the  central  lobe  of  the  platinum-centered  d 
orbital . 

The  origin  and  character  of  the  nodes  in  these  orbitals  was  clarified 
by  a  procedure  which  is  illustrated  In  Figure  7  for  (H^Pj^PtH^.  The  eigen¬ 
values  were  recalculated  for  new  geometries  in  which  bonds  of  interest 
were  selectively  stretched.  For  example,  increasing  the  Pt-H  distance  of 
L^PtHz  lowers  the  energy  of  the  6aj  orbital.  We  conclude,  then,  that  this 
orbital  Is  antibonding  with  respect  to  the  platinum-hydrogen  interaction. 

To  estimate  the  strength  of  the  bonding  (anti bond ing)  interaction  between 
two  nuclei,  we  relied  on  two  criteria.  One  was  the  magnitude  of  the  change 
in  energy  which  accompanied  a  small  change  in  bond  length  (Figure  7):  orbitals 
whose  energies  were  sensitive  to  these  changes  were  considered  relatively 
(anti ) bond ing.  The  second  was  the  amount  of  charge  indicated  to  lie 
between  the  nuclei  in  the  contour  diagrams  (Figures  2-4):  if  these  diagrams 
showed  a  significant  area  (volume)  at  the  highest  contour  level  (for  example, 
between  platinum  and  hydrogen  in  the  Vo^  orbital  of  l^Ptt^)  the  orbital  was 
considered  bonding. 

Having  identified  the  molecular  orbitals  responsible  for  bonding 
platinum  to  hydrogen  and  carbon  in  these  complexes,  it  is  possible  to 
interpret  orbital  correlation  diagrams  describing  the  conversions  L^PtXY  •+■ 
l^Pt  +  XY  (Figure  8).  These  calculations  contain  an  arbitrary  element:  as 
indicated  earlier,  we  do  not  know  the  LPtl  bond  angle  of  the  L^Pt  fragment 
produced  as  a  product.  Orbital  correlations  are  given  for  two  extremes: 


Figure  7.  The  response  of  the  eigenvalues  of  LjPtH^  (a)  to  a  10%  increase  in 
the  Z-coordinate  of  the  hydrogen  atom  (b).  These  shifts  in  energies  were  used 
to  classify  the  energies  as  bonding  (B) ,  non-bonding  (N)  or  antibondlnq  (A) 
with  respect  to  the  Pt-H  interactions. 
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one  having  the  L-Pt-L  angle  corresponding  to  the  starting  material 
(LLPtL  s  90°)  and  the  second  to  that  of  the  final,  relaxed,  product 
(LLPtl  -  180°).  The  information  summarized  in  Figure  6  indicates  that 

choosing  another  bond  angle  will  not  significantly  influence  the  conclusions, 
particularly  insofar  as  only  the  relative  rates  of  decomposition  of  the 
three  platinum  complexes  are  of  interest. 

Figures  2-4  and  the  orbital  correlation  diagrams  of  Figure  8 
contain  several  items  of  information  relevant  to  the  elimination  of  X-Y 
from  l^PtXY.  First,  the  orbital  correlation  diagrams  indicate  that  all 
three  eliminations  considered  are  allowed  by  orbital  symmetry.  Second, 
Figures  2-4  imply  that  the  strength  of  the  bonding  between  Pt  and  X(Y) 
is  not  due  to  a  single  orbital,  but  represents  the  sum  of  interactions 
in  several  orbitals.  Third,  these  figures  suggest  a  definite  difference 
in  the  distribution  of  Pt-X(Y)  bonding,  nonbonding,  and  antibonding 
contributions  among  these  orbitals.  For  L^PtCCH^)^  the  upper  orbitals 
have  Pt-CH^  bonding  character  and  there  are  no  occupied  Pt-CH^  anti  bonding 
orbitals.  For  l^PtH^.,  the  highest  occupied  orbital  is  Pt-H  antibonding, 
and  the  lower  ones  are  Pt-H  bonding.  For  l2Pt(H)CHj,  an  Intermediate 
situation  is  observed.  Thus,  these  Figures  contain  the  suggestion  that  the 
relative  rates  of  reductive  elimination  (that  is,  essentially,  the  relative 
thermal  instabilities)  from  these  complexes  correlate  with  the  occupancy 
of  Pt-X(Y)  antibonding  orbitals. 

L2Pd(CHj)2  and  L2Ni(CHj)2.  Figure  9  gives  eigenvalues ;  Figures  10 

^  <*V  A/  A/  /v  -W  *>+  /A*  A/  A/  A/  A/  A« 

and  11  summarize  information  concerning  wavefunct ions ;  Figure  12  gives 
orbital  correlation  diagrams  for  reductive  elimination.  All  are  given  in 
forms  analogous  to  those  used  with  l2Pt(CHj)2,  to  facilitate  comparison 


ENERGY  (eV) 


Is  In  its  ground-state  geometry. 
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The  Inferences  from  the  calculations  for  l^PdCCH^  and  L^Ni  (CHj)^ 
are  similar  tr  those  for  the  platinum  complexes.  The  orbital  correlation 
diagrams  indicate  that  all  of  the  reactions  are  allowed  by  orbital  symmetry. 
Examiration  of  the  occupied  orbitals  indicates  that  the  HOMO  of  L^Ni  (CH^ 
(the  8aj  orbital)  is  weakly  antibonding  with  respect  to  the  C-Ni  bond  (the 
orbital  is  predominantly  a  nickel  dz2  orbital,  with  relatively  little  elec¬ 
tron  density  on  carbon).  The  orbitals  of  (^PdCCH^Jj  and  LgPt^CH,)^  are 
qualitatively  indistinguishable:  neither  has  identifiable  occupied  M-C 
antibonding  orbitals. 

Discussion 

^ 

The  bonding  between  M  and  X(Y)  in  L^MXY  seems  to  be  distributed 
among  a  number  of  orbitals  of  similar  energies.  To  the  extent  that  the 
rates  of  reductive  elimination  from  these  complexes  are  determined  by  their 
electronic  structure,  these  rates  are  not  expected  to  be  entirely  controlled 
by  the  properties  of  the  HOMO,  or  by  orbital  symmetry  considerations,  but 
rather  by  the  collective  action  of  several  valence  orbitals.  Nonetheless, 
a  limited  number  of  these  orbitals  undoubtedly  are  particularly  important. 
Orbitals  with  the  character  of  1  or  2  are  clearly  important  bonding  and 
antibonding  orbitals;  molecular  orbitals  formed  from  combinations  of  d_^Z, 
d  2_  2,  s,  and  p  atomic  orbitals  on  the  metal  (3,4)  also  appear  from  the 
orbital  contour  diagrams  to  be  able  to  contribute  significantly  to  bonding 


or  antibonding. 
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One  clear  correlation  emerges  from  this  work  relating  molecular 

orbital  structure  to  rates  of  reductive  elimination.  Compounds  which  decompose 

rapidly  (l^PtH^,  l^NifCH^^)  have  occupied  orbitals  with  antibonding  M-X 

character;  compounds  which  decompose  slowly  (LjPtCCH^j)  have  occupied  M-X 

bonding  orbitals  and  vacant  M-X  antibonding  orbitals.  In  the  simplest  terms, 

this  correlation  is  almost  trivial:  it  Implies  that  compounds  in  which 

the  M-X(Y)  bonds  are  weak,  decompose  more  rapidly  than  those  in  which  these 
22 

bonds  are  strong  .  The  interesting  question  then  becomes:  what  is  the 
basis  in  atomic  or  molecular  structure  for  the  changes  in  molecular  electronic 
structure  which  are  reflected  in  the  occupancy  of  antibonding  orbitals? 


MiilfifijliiittaiHHH 


We  hypothesize  that  the  difference  between  the  energies  of  the  interacting 
orbitals  of  M  and  X(V)  Is  one  determinant.  Figure  13  summarizes  the  argument 
which  forms  the  basis  for  this  hypothesis  in  a  highly  simplified  and 


schematic  way.  Consider  a  group  of  orbitals  for  the  L^Pt  moiety  (|_  =  go°) 
of  l-2ptX2f  and  the  HOMO  (a_  symmetry)  for  the  X2  moiety.  Assume,  for 
simplicity,  that  the  orbital  interacts  only  with  the  I^Pt  of  the  correct 
symmetry  closest  to  it  in  energy.  (Strong  interaction  between  orbitals  of 
similar  energies  is  usually  rationalized  on  the  basis  of  orbital  size  and 
overlap'*'  .  In  this  irstance,  since  Pt  and  H  or  CH^  have  somewhat  different 
orbital  sizes,  this  assumption  is  not  strictly  justifiable.  Nonetheless, 
the  qualitative  argument  remains.)  If  it  is  an  occupied  orbital  (of  the 
correct  symmetry  for  overlap)  of  the  L2M  complex  which  matches  the  a x 


orbital  most  closely  in  energy,  the  resulting  M-X  bonding  and  antibonding 
orbitals  may  both  be  occupied  (Figure  13A,  exemplified  by  L2PtH2  an<* 
l^Nl  (CH^^)  •  If  it  Is  an  unoccupied  L^M  orbital  which  lies  closest  in 
energy  to  a ^  (Figure  138,  LjPtfCHj^,  LjPdfCHj^  or  if  the  HOMO  of  l^M 
interacts  strongly  with  a*  (Figure  13C;  note  the  6b_  level  in  L,Pt(CH,)_) 

^2  "4  ^  j  c 

the  cr£v  orbitals  resulting  from  these  interactions  should  be  vacant. 

M X2 

Thus,  in  brief,  if  the  orbital  electronegativities  are  similar  for  aY  and  the 

2 

highest-lying  I^M  orbital  of  the  same  symmetry,  the  MX2  fragment  can  be 
expected  to  have  an  occupied  anti  bonding  orbital  and  to  be  unstable  with 
respect  to  reductive  elimination;  if  these  orbital  electronegativities  are 
enough  different  that  one  of  the  strongly-interacting  pair  of  high-lying 
orbitals  is  vacant,  the  MX2  fragment  should  have  a  vacant  antibonding  orbital 
and  be  relatively  stable  to  reductive  elimination.  In  other  terms,  the 
interaction  suggested  by  Figure  13A  is  predominantly  covalent  while  that  in 


1  j 


LUMO 


Figure  13.  A  schematic  energy  level  diagram  summarizing  the  important 
interactions  in  l^MXY. 
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Figures  138  and  13C  has  more  or  less  ionic  character. 

Reference  to  Figures  2,  4,  10  and  11  indicate  that  the  important 
elements  in  these  orbital  correlation  diagrams  are  approximated  by  the 
interactions  outlined  in  Figure  13,  although  they  are  obviously  more 

I 

complicated,  since  several  L^M  orbitals  have  the  correct  symmetry  to  interact  ; 

with  aY  and  a*  . 

2  *2 

Inspection  of  calculated  charges  in  the  several  complexes  examined  j 

here  supports  these  hypotheses  concerning  the  importance  of  electronega¬ 
tivities  (Table  1).  For  the  series  the  complex  in  which  the 

least  charge  transfer  from  carbon  to  metal  would  take  place  during  reductive 
elimination  is  L^NitCHj)^;  this  complex  is  also  the  one  for  which  reductive 
elimination  is  most  rapid.  For  the  series  I^PtXY,  the  complex  for  which 

the  least  charge  is  transferred  to  the  metal  during  reductive  elimination  i 

is  LjPtH^.  This  complex  is  also  the  one  which  reductively  eliminates 
most  rapidly  of  the  series  containing  platinum 

This  model  correlates  the  rates  of  reductive  elimination  with 
the  presence  of  occupied  M-X(Y)  antibonding  orbitals,  and  by  inference  with 
the  relative  orbital  electronegativities  of  the  orbitals  and  the 
"stretched"  a^y  orbital.  It  hypothesizes  that  the  relative  electronic 
energies  of  the  starting  materials  determine  the  rates  of  these  reductive 
elimination  reactions.  A  number  of  other  possible  models  for  these  reactions 
might,  of  course,  also  be  considered.  Several  of  these  follow.  We  note  that 
each  of  the  factors  emphasized  in  these  models  undoubtedly  contributes 

to  some  extent  to  the  rates  of  these  reactions  and  that  they  are  not  j 

| 

necessarily  orthogonal. 
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Table  1  Calculated  Atomic  Charges  In  Complexes  of  the  Structure  LjPtXY 

(L  -  PH3,  X,  Y  -  H,  CH3) 


Complex 

P 

M 

C 

H 

4PtH2 

+0.28 

-0.20 

-0.04 

L2PtHCH3 

+0.298 

-0.160 

-0.107 

-0.024 

L2Pt(CH3)2 

+0.28 

-0.001 

-0.176 

L2Pd(CH3)2 

+0.27 

-0.061 

-0.140 

L2Ni(CH3)2 

+0.44 

-0.59 

-0.013 

H2 

0.0 

(Hc),C-H 

+0.174 

-0. 129— 

W 

+0.177 

-0.067 

L2Pt(90°) 

+0-393 

-0.439 

L2Pt(l02°) 

+0.403 

-0.460 

L2Pt(I80°) 

+0.422 

-0.521 

L2Pd(102°) 

+0.525 

-0.559 

L2Ni(l02°) 

+0.372 

-0.601 

a  ° 

—The  value  for  Is  +0 .013-  The  C-H  distance  assumed  is  2.73  A,  and  corresponds 

to  the  stretched  value  in  L2Pt(H)CH3.  —The  assumed  C-C  distance  is  2.72 

and  corresponds  to  the  value  in  L,Pt(CH3)2. 
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1)  Exothermic! ty.  The  react (on  rate  for  the  decomposition  of  I^MXY 

might  correlate  with  the  free  energies  of  these  reactions:  such  correlations 

are  numerous  in  kinetics,  although  they  usually  apply  only  to  a  restricted 

27 

range  of  variations  in  the  reactant  .  The  arguments  presented  in  this  paper 

deal  with  only  one  contribution  to  this  free  energy:  viz ,  the  enthalpy  of  the 

reactants.  The  enthalpies  of  the  products  are  the  second  contributor,  and 

we  observe  that  the  bond  energies  of  the  products  (kcal  mol  ^ “111  > 

“  99  >  *  83)  also  correlate  with  the  rates  of  decomposition  of 

L^PtXY.  The  relative  rates  of  decomposition  might,  in  principle,  be  dominated 

by  the  relative  energies  of  either  the  starting  materials  or  the  products. 

The  data  needed  to  identify  the  more  important  contributor  are  not  available 

in  the  series  examined  here.  Nonetheless,  comparison  of  our  results  with 

9 

those  obtained  by  Hoffmann  et_  a]_?  using  extended  Huckle  methods  suggests 
that  in  at  least  some  areas  the  response  of  the  transition  state  to  changes 
in  structure  may  be  more  like  that  of  products  than  reactants  (see  below). 

2)  Relief  of  Steric  Strain.  Changes  in  nonbonded  interactions  within 

crowded  transition  metal  complexes  are  known  to  be  of  major  importance  in 

determining  the  rate  of  many  organometal 1 ic  reactions,  including  reductive 

2  8 

el iminat ions  and  oxidative  additions  ’  .  The  differences  in  size  of  H  and 
CHj,  and  of  Pt(ll),  Pd(ll),  and  Ni(ll),  could  in  principle  result  in  contri¬ 
butions  to  rates  from  nonbonded  interactions  between  X(Y)  and  L  larger  than 
those  due  to  differences  in  bond  energies. 


Of  the  many  such  effects  which  could  be  imagined  a  priori,  one  in  particular 
is  suggested  by  the  contour  diagrams  of  Figures  2-4.  During  reductive  elimi¬ 
nation,  the  X-Y  group  forms  a  bond.  The  l s_ orbital  on  hydrogen  is  spherically 
symmetrical;  the  £  (or  sg^)  orbitals  on  carbon  are  not.  Thus,  three-center 
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bond  I  ng  of  the  type  5  might  be  more  important  for  hydrogen  (6)  than  methyl  (7) 

A/  ^  <V 

The  6aj  of  l2PtH2  shows  sc  -''idence  of  H-H  bonding 


5  6  7 


in  the  ground  state;  a  corresponding  concentration  of  charge  between  the 
carbon  atoms  of  l^PtCCH^,  or  between  the  carbon  and  hydrogen  of  L2Pt(H)CHj, 
is  not  evident. 

9 

Brief  comparison  of  our  results  with  those  of  Hoffmann,  et_  aK  is  useful 

Hoffmann  reached  three  major  conclusions:  (1)  The  more  strongly  electron 

donating  is  X  (and  Y) ,  the  more  rapid  is  reductive  elimination;  (2)  The  more 

strongly  electron  donating  is  L,  the  slower  is  reductive  elimination;  (3)  The 

barrier  to  reductive  elimination  Is  determined  In  major  part  by  b2  orbitals  of 

type  i.  The  first  two  of  these  conclusions  have  some  experimental  support. 

2  8 

In  general,  the  complexes  LjMXj  ar®  more  stable  for  X  -  CF^  ,  Cl,  Ph  than 
for  X  •  CHj ;  Halpem,  et_  ah  have  reported  that  the  rate  of  loss  of  methane 
from  (X-Ph^P)2Pt (CHj)H  is  more  rapid  when  X  Is  electron-w? thdrawing.  Our 
calculations  apply  only  to  the  electronic  structures  of  the  reactants. 


A  large  (several  eV)  lowering  in  the  energy  of  the  X(Y)  orbitals,  or  a  large 
Increase  in  the  energy  of  the  orbitals  should  increase  the  probability 
that  PtXY  antibonding  orbitals  are  occupied,  and  predict  a  destabilization 
of  the  complexes.  The  influence  of  smal 1  changes  of  the  type  encountered  in 
substituents  is  not  obvious.  Qualitatively,  since  charge  density  on  the  metal 
center  of  L^M  is  larger  than  that  in  LjMXY,  electron  donating  L  should  de¬ 
stabilize  the  product  more  than  the  reactant.  Similarly,  since  reductive 
elimination  involves  charge  transfer  from  X (Y)  to  metal,  the  more  electron- 
withdrawing  X(Y),  the  more  stable  should  be  l^MXY.  Thus,  the  results  of  our 
studies  provide  an  adequate  ad  hoc  rationalization  of  observed  substituent 
effects,  and  can  be  interpreted  to  be  consistent  with  the  potential  surfaces 

calculated  by  Hoffmann,  et  al_.  They  do  not,  however,  lend  themselves  easily 
to  predictions  of  rate  differences  resulting  from  small  perturbations  in 
substituent  structure.  The  third  conclusion  of  Hoffmann,  that  the  energy  of 
an  orbital  of  type  1  dominates  the  changes  in  energy  in  going  from  reactant 
to  transition  state,  cannot  be  compared  with  our  results,  since  we  do  not 
calculate  total  energies. 

Thl’s  work  was  suPPorted  by  NSF  (OMR -78-241 85)  anij  the 
Office  of  Naval  Research. 
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Captions 

Figure  1.  Eigenvalues  for  l^PtXY  [XY  -  HH  (top),  CH^H  (middle),  CHj-CHj 
(bottom)].  Only  the  orbitals  of  L2Pt  and  XY  making  major  contributions  to 
the  moiecuiar  orbitals  of  l^PtXY  are  indicated.  The  eigenvalues  for  XY  are 

O 

those  calculated  at  the  distances  assumed  in  the  complex  (r^  *  2.80  A, 

o  a 

rCH  -CH  "  2.73  A,  r^  ■  2.72  A).  In  addition,  the  eigenvalues  for  CH^-CH^  are  also 

calculated  using  the  bent  geometry  (C2y  symmetry)  this  moiety  displays  in  the 
L2Pt(CH^)2  complex.  Calculations  of  the  distorted  CH^-H  molecule  displaying 
both  Cjy  and  C^  symmetries  revealed  insignificant  differences.  The  figure 
displays  the  eigenvalues  obtained  from  the  CH^-H  fragment  having  C^v  symmetry.) 

I  -  PH^  in  each  case. 

Figure  2,  Orbitals  important  to  bondinq  in  L2PtH2-  The  first  column  represents 
the  simplified  interpretation  of  the  bonding  in  the  orbital.  The  representation 

is:  . -  - —  stronaly  bonding; - weakly  bonding; 

strongly  antibondinq.  Thus,  for  example,  the  2b2  orbital  has  strong  Pt-H 
bonding  character  and  weak  Pt-P  bonding.  The  two  hydroqens  (and  the  two 
phosphorus  atoms)  are  weakly  antibonding  with  respect  to  one  another.  No  line 
between  nearest  neighbor  atoms  (as  between  Pt  and  P  in  the  4a^  orbital)  indicates 
a  weak  or  non-bonding  interaction.  If  a  center  is  omitted  from  the  diagram,  the 
charge  assigned  to  it  is  considered  sufficiently  small  that  it  effectively  does 
not  participate  in  the  orbital.  The  second  column  sketches  the  orbital  makeup 
in  terms  of  atomic  orbitals.  The  numbers  represent  the  fractional  charge 
(expressed  as  a  percentage  of  two  electrons)  centered  on  each  center.  These 
charges  do  not  add  to  100  because  the  charge  on  the  hydrogens  in  PH^  and  CH^ 


have  been  omitted  from  the  diagram.  The  major  contributions  of  the  orbitals 
on  platinum  are  indicated  on  the  diagram.  For  example,  the  2a^  orbital 
contains  important  contributions  from  d  ^  2  anc*  —  orbitals  centered  on 

platinum.  The  last  column  contains  a  section  through  the  wave  function.  The 
contours  are  plotted  at  values  of  *0.003,  *0.009,  *0.027  and  *0.081,  starting 
from  the  outermost. 

Fiaure  3-  Orbitals  important  to  bonding  in  L^PtCHjCH^.  The  caption  for 
Figure  2  contains  nomenclature. 

Figure  k.  Orbitals  important  to  bonding  in  LjPtfCH^jj-  The  caption  for 
Figure  2  contains  nomenclature. 

Figure  5-  Orbitals  for  L^Pt (0)  for  angle  LPtL  ■  141°.  The  caption  for 
Figure  2  contains  nomenclature. 

Figure  6.  Eigenvalues  for  L^PtfO)  as  a  function  of  anqle  LPtL. 

Figure  7-  The  response  of  the  eigenvalues  of  LjPtH^  (a)  to  a  10$  increase  in 
the  Z-coordinate  of  the  hydrogen  atom  (b).  These  shifts  in  enerqies  were  used 
to  classify  the  energies  as  bonding  (B),  non-bonding  (N)  or  antibondinq  (A) 
with  respect  to  the  Pt-H  interactions. 

Figure  8.  Orbital  correlations  diagrams  for  decompositions  of  L£ptXY  to 
XY  and  L^Pt  (anqle  ■  102°  and  l80°).  XY  is  in  its  ground-state  geometry. 


Figure  9.  Eigenvalues  for  [M  -  Pt,  Pd,  NI].  Only  the  orbitals 

of  L2M  making  major  contributions  to  the  molecular  orbitals  of  L2M(CHj)2 
are  indicated.  L  ■  PH^  in  each  instance. 

Figure  10.  Orbitals  important  to  bonding  in  I^NiCCH^.  The  caption 
for  Figure  2  contains  nomenclature. 

Figure  11.  Orbitals  important  to  bonding  in  I^PdfCH^.  The  caption  for 
Figure  2  contains  nomenclature. 

Figure  12.  Orbital  correlation  diagrams  for  decomposition  of  L2M(CH^)2  to 
C2H^  and  LjM  (angle  ■  102°).  C2Hg  Is  in  its  ground-state  geometry. 

Figure  13-  A  schematic  energy  level  diagram  summarizing  the  important 
interactions  in  L2MXY. 


